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Abstract: The "glassy" superconducting transition at high magnetic fields can be induced 
by columnar disorder. A model is proposed in which the thermodynamics of Bose condensa- 
tion of Cooper pairs into the lowest Landau level eigenstate of the random potential can be 
solved exactly. The solution reflects a peculiar character of the high-fleld limit: For example, 
the effective dimensionality of the transition is shown to be a function of magnetic fleld. 
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The problem of superconducting transition in the presence of strong disorder is of both 
practical and theoretical interest. Technologically, the goal is to introduce defects into the 
sample in a way that maximizes pinning of vortices and increases critical currents 0]. The 
theoretical challenge is to understand the mechanisms and the nature of superconducting 
transition for various types of disorder. A variety of novel phases have been proposed, dif- 
fering in the cases of point like 0, ^ and line-like disorder 0, ^ . In this Letter we present a 
theory of superconducting transition at high magnetic fields (> 1 Tesla in HTS) induced by 
the presence of columnar (line-like) defects. In the absence of disorder, the high-field fluc- 
tuations of the order parameter, ip{r), are strongly enhanced by formation of Landau levels 
(LLs) for Cooper pairs Such fluctuations lead to D D — 2-dimensional reduction in 
the pairing susceptibility, Xsc{r',r'), and eliminate the superconducting (Abrikosov) transi- 
tion for D = 2,3 0]. The Abrikosov phase is then replaced by a new fluctuation-induced 
state, the density-wave of Cooper pairs (SCDW), in which the thermal average {\'ip{r)\'^) 
has a weak modulation accompanied by only a short range phase coherence 0, P|. In the 
presence of disorder the LL degeneracy is lifted and a possibility of superconducting tran- 
sition is restored. Xsc can now diverge at some finite temperature, Tsc{H), determined by 
the strength of disorder 0. At T^dH) the normal state is unstable to Bose condensation 
of Cooper pairs into the lowest energy eigenstate of the random potential, which we argue 
extends over the whole sample in situations of experimental interest. Furthermore, for ex- 
perimentally relevant parameters, Tsc{H) can be far above the SCDW transition line over 
much of the H — T phase diagram, allowing us to treat the correlations that produce SCDW 
in an approximate way. 

We consider a realistic model for a superconductor in a magnetic field parallel to the 
columns which is exactly solvable. The model exhibits "dimensional transmutation", i.e. 
the effective dimensionality of the transition changes continuously as a function of magnetic 
field. This effect is a direct consequence of analytic properties of LL wavefunctions and is a 
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signature of the high- field hmit. We determine the transition hne m H — T phase diagram, 
the Edwards- Anderson order parameter, and the behavior of correlation length, specific heat 
and magnetic susceptibility in the vicinity of the transition. There are similarities between 
the transition considered here and the one in the spherical model for spin-glasses 

We are interested in strongly anisotropic layered superconductors described by the 
Ginzburg-Landau (GL) Lawrence-Doniach model, with magnetic field perpendicular to the 
layers. Fluctuations of the magnetic field are neglected (k ^ 1). We focus on the high-field 
limit, where the LL structure of Cooper pairs dominates the fluctuation spectrum: This is 
the case for fields above Hi, ^ (6'/16)//c2(0)(T/Tc(0)), where 6 is the Ginzburg fluctuation 
parameter IlO|. (For instance, in BSCCO 2:2:1:2, 9 ^ 0.045 and ^ I Tesla.) In this 
regime, the essential features of the physics are captured by retaining only the lowest Lan- 
dau level (LLL) modes. This is the renormalized GL-LLL theory [|^, |T0[. The partition 
function is Z = J D[iIj*iIj] exp(— 5), and 
d ^'^ r ( 



n=l ■ 



where a'(T, if) = a(T — Tc2{H)), d is the effective layer separation, n is the layer index, A > 
is the effective strength of the defects and a, f3 and rj are phenomenological parameters. The 
magnetic field is assumed to be parallel to columnar defects, the effective potentials of all 
defects the same and well-represented by delta-functions. Random variables in the problem 
are 2D coordinates of defects, {rl}. We assume that columns of damaged superconductor are 
distributed according to the Poisson distribution Pj^{fl, ...rj^) = {e^P^p^)/N\ where is 
the probability for finding N impurities at the positions fl, ...rjy, A is the area of the system 



and p is the concentration of impurities [11]. After rescaling the fields and the lengths as 



{2d/32TrP /ry/^ip i), r/(/v^) ^ r, where / is the magnetic length for charge 2e, the 
quartic term can be rewritten 

-j: J d'r\^n\' = ^j:pA{n){j d'fi^r^ir (2) 



where N = A/2'Kf is the degeneracy of the ILL and /5a (^) = / |'^nH/(/ I'^nP)^ is the 
generahzed Abrikosov ratio corresponding to configuration ipni'f^)- We now observe that 
/3a (^) is only weakly dependent on the actual configuration, the well known example being 
the small difference in /?a between triangular and square lattice of zeroes llO|. Thus, we may 
substitute PAi^) in the quartic term by its thermal average, (Pa), and treat this as an input 
to the theory. This approximation neglects the non-perturbative lateral correlations that 
produce the SCDW transition [0. It is justified if TgdH) is far above the SCDW transition 
line. In that case the SCDW correlations enter only very close to the transition and can be 
ignored in most realistic situations. Since superconducting and SCDW transitions arise from 
two distinct mechanisms, the respective transition lines scale differently in the H — T phase 
diagram and, for moderate disorder, we are assured of a wide region near Hc2{T) where the 
neglect of SCDW correlations should be justified (see Fig. 1) [|T^. 

After the /5a (^) ^ (Pa) substitution the thermodynamics of the model becomes exactly 
solvable. We first introduce variables to decouple the quartic term and integrate over 
the fields {ip*:i^)- This leads to Z = / Hn '^^n exp(— A^S"), where 

n \Pa) Jo 

(3) 

We drop terms coming from the rescaling of V'n(^) and introduce dimensionless combinations 
of GL parameters gr^^a,x = {ri,a', \/2ttP} x ^ (dnP)/ (T/3). The density of states for a 
disordered system in the LLL can be found exactly by using the supersymmetric formalism 
13| . For the Poisson short-range scatterers it is given by 



pJV) = -Im-^ln r dtexp(m - f f —(1 - e"'^)) (4) 
TT dV Jo Jo y 

where / = p27r/^. In the thermodynamic limit N ^ oo, the partition function is completely 
determined by the saddle-point of S'. Assuming that the saddle point is at x independent 
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of the layer index, we finally write the free energy above the critical temperature 
F x2 



^J_^dkJ^ dVpf{V)\n[gr,e{k)+g^ + x + gxV] , (5) 



NNlT {Pa) 2 
where e{k) = 1 — cos(A;) and x is determined by the solution of 

4 J-i Jo gne{k) + ga + X + gxV 
In Eq. 6 it is important to know the behavior of density of states at low energies. 
For / < 1, density of states has a delta-function singularity at = 0, whilst for / > 1, 
Pfiy) ~ V-^~'^ when V ^ [T^. When V < 0, p{V) = 0, as also can be inferred from Eq. 4 
jTSl . The transition line, TgdH), in the H—T diagram is determined by Eq. 6 and x+ga = 0, 
which corresponds to condensation of Cooper pairs into k = and = eigenstate of the 
random potential. It is easily seen that there will be a non-zero transition temperature only 
if concentration of impurities and magnetic field are such that / > 3/2. Below this value of 
/ LLL degeneracy is not sufficiently lifted by the random potential and thermal fluctuations 
prevent a finite temperature phase transition in our model [Q. Experimentally, this should 
manifest itself as a drop in transition temperature when the field exceeds a certain value. 
/ = 3/2 determines the effective lower critical dimension for our model. After introducing 
dimensionless quantities t = T/Tc(0), h = H/Hc2{0) and A' = XHc2{0)/(j)oaTQ, where 0o is 
the flux quantum, we perform the integration over wave-vector k in Eq. 6 to obtain the 
expression for transition temperature 



tscih) = il-h) 



2A' Jo ./y2 + {2t^V) /{hX'aT^ 



(7) 



Notice that when A' ^ we have tgcih) — * 0, while for A' oo, tg^ih) ^ 1 — /i [|15 
Also, with increasing parameter /, tsc{h) increases. This is related to the observation in Ref. 
1 that the irreversibility line shifts to higher temperatures with increasing doses of irradiation 
with heavy ions. Numerical solution for tsc{h) is displayed in Fig. 1 for A' = 1, rj/aTo = .01, 
9 = 0.03, (Pa) = 1-3 and / = 0.04/ h. We have set {Pa) to a constant for simplicity. If 



Hc2{0) ~ 100 Tesla, / = 0.04/ h corresponds to average distance between defects of 225A (at 
1 Tesla, I = 180A). The diameter of the columns depends on the size and energy of particles 
used for irradiation but it is about 50A and hence much smaller than the magnetic length 
for the fields of interest. Thus, representing defects by delta-functions is appropriate. 

As temperature drops below tsc{h), x remains at the value it had at the transition. 
There is now macroscopic occupancy of the lowest energy state at ^ = and k = 0. As 
is well known, condensation into this state is possible only if the state is extended. It is 
a special feature of this problem that this indeed is the case for certain range of impurity 
concentration. The density of states, Eq. 4, changes at \^ = from being infinite when 
/ < 2, to being zero when / > 2. Thus, for fields and impurity concentrations such that 
parameter / < 2, true extended states (which always exist in the LLL [T^) must lie at 
the bottom of the impurity band, since the number of states there diverges. The change of 
behavior in the density of states at / = 2 could be caused by the fact that the mobility edge 
shifts to positive energies at some /o > 2, leaving spread-out but localized states at \^ = 0, 
which now becomes the tail of the distribution. Numerical diagonalization studies indicate 



that mobility edge is indeed located near the band center for / > 4 |[T^. Thus, strictly 
speaking, our model is appropriate for / < /q. However, even for / above but close to /o, 
which is often the case for fields and concentrations of experimental interest, the states at 
V = are still near mobility edge and will appear extended in a finite size sample. On this 
basis, we expect that useful information about the transition can still be obtained within 
our model. 

With these cautionary remarks in mind, the natural order parameter is the thermal 
average of the component of ipni''^) corresponding to the eigenvalue with V = and k = 0. 
This is {ipo,o) = {NNL{ga\t=tscih) — da)/ {Pa))^^^- The disorder average value of the field is 



('0„(r)) = J2v,k 6xp (ikn) (V'y.fc) = 0, due to random phases of the state 0v'=o(^- Under 



the assumption that the lowest state is extended through the sample, |0y=o(r)p ~ 1/-^; the 
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Edwards- Anderson order parameter |T8| Qea = then equals 



lEA - J-TT i9a\t=tscih) - 9a) (8) 

\Pa) 

below tsc{H), and is zero above. Thus, qea oc itsdh) — t)^^, with the exponent /? = 1/2. 
The free energy below tsdh) is 

To calculate the exponents that determine the divergence of correlation lengths parallel 
and perpendicular to the field we first note that from Eq. 6 and the definition of the critical 
line it follows 



{9a + X) 



^+ 4 7-/^0 {9Ak)+9xV){g,e{k)+g^^ + g^ + x)\-^'' 9a\t=u.iH). 

(10) 

The integral in the last equation diverges for / < 5/2 as {ga + x)^^^^"^ when the transition 
line is approached from above, and it is finite for / > 5/2. Thus, we obtain {ga + x) oc 
{t — tscY^^^"^^'^'' for /< 5/2 and {ga + x) oc (t — tgc) for /> 5/2. The same behavior follows if 
the transition line is approached along the line of constant temperature. This determines the 
value of the exponent z/y = 1/(2/— 3) for / < 5/2 and the classical value = 1/2 for / > 5/2, 
where the correlation length parallel to the field is oc [t — tsc(^)]~'^"- The concentration 
corresponding to / = 5/2 determines the effective upper critical dimension in the problem. 
We now turn to the correlation length perpendicular to the field, ^± oc [t — tsc{h)]~'^-^ and 
study the susceptibility associated with Edwards- Anderson order parameter. 



Xi?vi(r - r') = (V^*(r)?/'„(r'))(V'„(f)V^*(r')) . (11) 

After expanding the field operators in the eigenbasis of random potential we obtain 

_ r') = / dV^dV,dk^dk,F{r-r>,V,,V,) 

J (g^ki) + gxVi + g^ + x){gAk2) + gxV2 + ga + x) 
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where the function F is the two-particle spectral density |]T9| 



F{r- r', Ki, V2) = E ^(^1 - yi)^{y2 - l^,)0*(r)0,(r')0,(r-')0*(rO (13) 

and (/)j(r) are the eigenstates of the random potential. If we now introduce V = [Vi + V2)/2 
and uj = {Vi — V2) /2, for V close to the mobility edge and small {q, u), the Fourier transform 



of F has a diffusive form iT9l EOl 



F{q,V,uj) = in2f 2/ \\ ' 

7r(t<;'' + q^D'^[q^/uj)) 

where D{q^ /uj) is the generalized "diffusion constant". Assuming this form for F{q, V, u) and 
rescaling everything by the appropriate power of temperature in Eq. 12, we obtain z/_|_ = z/y. 
Note that the density of states PfiV) is roughly constant except in a narrow region, typically 
1% of total bandwidth, around V = 0, where it either diverges or vanishes. Hence, unless 
one experimentally probes the system very close to the transition, the observed correlation 
length exponent would be the one corresponding to f = 2, i.e. = v\\ = 1. This agrees 
well with experimental results of Ref. 5. 
Magnetization per unit volume equals 

M 2T^Vth ( 2x \ 



Below the transition line this coincides with the usual mean-field result. Above the transition 
line qEA = and from Eq. 10 it follows that at constant temperature close to the transition 
(da + x) (X [h — hscit)Y^^^^^^'^^ when / < 5/2 and (ga + x) oc [h — hsdt)] otherwise. Thus the 
magnetic susceptibility is a smooth function of the field at the transition for / < 2, but has an 
upward cusp for 2 < / < 5/2 and a discontinuity for / > 5/2. The size of this discontinuity 
depends on the location of the transition in the H — T diagram. Differentiating the free 
energy twice with respect to temperature one obtains the specific heat. It is straightforward 
to show that at the transition it behaves the same way as susceptibility; smooth for / < 2, 
has a cusp for 2 < / < 5/2 and has the usual discontinuity for /> 5/2. More precisely, both 
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magnetic susceptibility and specific heat behave as [t — tsdh)]^" for 3/2 < / < 5/2, where 
a = {f — 5/2)/ (/ — 3/2). The behavior of the specific heat, order parameter and correlation 
length in our model is analogous to the one obtained from 0(2A^) vector model in the limit 
^ oo and in the effective dimension -De// = 2/ — 1 The magnetic susceptibility 

however, behaves differently at the transition; whilst it diverges in the 0(2A^) vector model 
with the exponent 7 = [{Deff/2) — 1]^^, it is finite in our case even below the effective upper 
critical dimension. This is a consequence of a diamagnetic nature of magnetization in our 
problem. 

We should stress again that the critical behavior of our model does not describe "true" 
critical properties of the GL-LLL theory with disorder, since we have ignored lateral SCDW 
correlations. Such correlations will always become important sufficiently close to the transi- 
tion. However, as it is clear from Fig. 1, there is a wide region in the H — T phase diagram 
where the superconducting transition lies far above the SCDW transition line for clean sys- 
tems. In this region, the "true" critical behavior will set in only very near the Tsc{H) line 
and our model should be appropriate in most experimental situations. 

In summary, we have studied the high-field superconducting glassy transition induced 
by columnar disorder. Using the exact density of states for random array of short-range 
scatterers in the LLL level and the assumption that the lowest eigenstate of such a potential 
is extended over a finite size sample under certain conditions, we have obtained the Edwards- 
Anderson order parameter, correlation length, magnetization and the specific heat close to 
the transition. The transition line in if — T phase diagram has also been calculated. The 
critical exponents are found to depend on magnetic field. 

This work has been supported in part by the David and Lucile Packard Foundation. 
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Caption: 

Figure 1: The H — T phase diagram for strongly type-II superconductor with columnar 
disorder {h = H/Hc2{0), t = T/Tcq). The full line represents the second-order phase transi- 
tion boundary between normal and "glassy" superconducting state for the set of parameters 
given in the text. Dashed-dotted line is the SCDW transition in clean system. Dashed hue 
is the mean-field hc2{t)- The LL approximation breaks down in the shaded region at the 
bottom. 
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